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Membrane ultrafiltration (UF) is a pressure driven process allowing for the separation and enrich¬ 
ment of protein solutions and dispersions of nanosized microgel particles. The permeate flux and the 
near-membrane concentration-polarization (CP) layer in this process is determined by advective- 
diffusive dispersion transport and the interplay of applied and osmotic transmembrane pressure 
contributions. The UF performance is thus strongly dependent on the membrane properties, the hy¬ 
drodynamic structure of the Brownian particles, their direct and hydrodynamic interactions, and the 
boundary conditions. We present a macroscopic description of cross-flow UF of non-ionic microgels 
modeled as solvent-permeable spheres. Our filtration model involves recently derived semi-analytic 
expressions for the concentration-dependent collective diffusion coefficient and viscosity of perme¬ 
able particle dispersions [Riest et al, Soft Matter, 2015, 11, 2821]. These expressions have been 
well tested against computer simulation and experimental results. We analyze the CP layer prop¬ 
erties and the permeate flux at different operating conditions and discuss various filtration process 
efficiency and cost indicators. Our results show that the proper specification of the concentration- 
dependent transport coefficients is important for reliable filtration process predictions. We also 
show that the solvent permeability of microgels is an essential ingredient to the UF modeling. The 
particle permeability lowers the particle concentration at the membrane surface, thus increasing the 
permeate flux. 


I. INTRODUCTION 

Microgels are colloidal dispersions of nano- to micro¬ 
sized visco-elastic gel particles (typically cross-linked 
polymer networks) that can swell or deswell in a con¬ 
trolled way in response to external stimuli such as tem¬ 
perature, pH value, salinity and particle concentration 
m- The microgel phase behavior is richer than that 
of classical colloidal hard spheres due to their size vari¬ 
ability and viscoelasticity. This makes microgels an ideal 
model system to study, for example, the effect of intrinsic 
particle viscoelasticity on fundamental physical processes 
like crystallization kinetics and glass formation |4]. Mi¬ 
crogels are used in a wide range of applications including 
drug delivery and advanced oil recovery, and for pharma¬ 
ceutical or cosmetic products 0H5]. 

In the production process of microgels, the reactor ef¬ 
fluent consists of a complex mixture of microgels and 
unreacted polymers dispersed in a solvent. Thus, various 
separation and concentration steps are necessary. Mem¬ 
brane UF is a suitable method to remove traces of free 
polymers and to concentrate the sample for later prac¬ 
tical applications [1]. UF is a pressure-driven process 
widely used in commercial large-scale plants by which 
smaller Brownian particles (like nanosized colloids and 
proteins) are concentrated and separated from smaller 
ingredients El 13- Its major advantage is the low energy 
requirement. UF has been traditionally applied to waste 
water treatment [HE], fruit juice concentration and clar¬ 
ification E3. and protein concentration nmg. Inside 
the human body, UF takes place in Bowman’s capsule at 
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the kidneys, where water and other small molecules are 
separated from blood m- 

UF must be distinguished from membrane microfiltra¬ 
tion of larger micro-sized particles which are also pres¬ 
sure driven. However, the applied transmenrbrane pres¬ 
sures are smaller and the permeation fluxes are sub¬ 
stantially larger than in UF. Moreover, and different 
from UF, Brownian motion and transmembrane osmotic 
pressure are small effects only. Consequently, shear- 
induced particle diffusion and migration are the leading 
hydrodynamic mechanisms in microfiltration. For cross- 
flow setup in microchannels, shear-induced migration en¬ 
hances the back migration away from the membrane, thus 
reducing the formation of a stagnant cake layer HMDS- 
Brownian particle motion dominates in UF, in contrast to 
microfiltration, with the local osmotic equilibrium being 
only weakly perturbed by the pressure driven flow. 

The efficiency of UF is influenced by the CP layer, 
which is the particles-enriched region next to the mem¬ 
brane surface created by advective-diffusive transport 
mna. The CP layer formed by mobile particles effec¬ 
tively increases the resistance to solvent transport and 
the local osmotic pressure, hence reducing the permeate 
flux. 

Membrane fouling is an undesirable irreversible mod¬ 
ification of the membrane caused by specific physical 
and chemical interactions between the membrane and 
the particles Q2iH22]. It gives rise to a stagnant cake 
layer deposited on the membrane and to the intrusion or 
clogging of particles inside the membrane pores. Foul¬ 
ing substantially lowers the permeate flux because of the 
added hydrodynamic resistance and can also change the 
membrane selectivity. Moreover, the filtration channel 
becomes constricted by the growth of the cake layer. In 
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order to avoid or reduce fouling effects in large volume 
concentration processes, UF is often performed in a cross- 
flow setup where, under steady-state conditions, the feed 
dispersion flow is directed parallel to the membrane with 
inlet and outlet ports. 

For small values of the applied transmembrane pres¬ 
sure (TMP), the permeate flux is found to be propor¬ 
tional to the TMP, while at larger TMP the flux increases 
sublinearly. Eventually, a limiting flux is reached, in¬ 
sensitive to any further increase in the TMP. Blatt et 
al. investigated the permeate flux in cross-flow UF us¬ 
ing thin-film theory (gel-layer model) under conditions 
of limiting flux [6i. The thin-film theory has been subse¬ 
quently improved within a boundary layer approximation 
by regarding explicitly the CP layer mechanism S5Vl25| . 
In these later works, the concentration dependence of 
the collective diffusion coefficient and/or the dispersion 
viscosity was considered using phenomenological expres¬ 
sions for the transport coefficients, with applications to 
specific systems like BSA proteins, sucrose, and dextran. 
Another model refinement was obtained by including the 
osmotic pressure buildup at the membrane surface op¬ 
posing the influence of the applied TMP mi na iizi- 

in addition, numerical solutions of coupled advection- 
diffusion, Navier-Stokes, and continuity equations by fi¬ 
nite difference methods have been used to analyze CP in 
cross-flow UF and hence to predict the permeate flux and 
the particle concentration profile at the membrane sur¬ 
face f?Sf!iI?l . In most of these works, the concentration- 
dependent osmotic pressure, collective diffusion coeffi¬ 
cient and dispersion viscosity were taken into account 
using approximate expressions for the special cases of 
solvent-impermeable colloidal hard spheres and charged- 
particles models describing sucrose, polyvinyl alcohols, 
and BSA and lactoferrin proteins. Typically employed 
approximate expressions for the sedimentation coeffi¬ 
cient have been the truncated virial expansion result by 
Batchelor [33] and the cell model result by Happel [34] . 
The semi-phenomenological expressions by Eilers and 
Krieger-Dougherty US have been used to describe the 
zero-frequency viscosity of neutral solvent-impermeable 
hard-sphere dispersions. The inaccuracy of these trans¬ 
port coefficient expressions has its impact in the respec¬ 
tive UF model predictions. 

Recently, Bouchoux et al. studied dead-end filtration 
of milk casein micelles [55] , These soft, deformable and 
solvent-permeable micelles can be regarded to some ex¬ 
tent as biological siblings of the synthetic microgels. One 
of the major issues of their work was to determine the 
permeability of the dispersion, which was obtained by 
a combination of osmotic stress and filtration measure¬ 
ments. Different from cross-flow, the dead-end or frontal 
filtration is an intrinsically non-steady process where the 
thickness of the immobile cake layer increases during the 
process. 

The present work includes the first comprehensive the¬ 
oretical study of particle permeability effects in the con¬ 
centration process of non-ionic microgel dispersions by 


cross-flow UF. The microgel particles are modeled as 
solvent-permeable spheres interacting by a hard-core po¬ 
tential. In contrast to previous works, we use accurate 
analytic expressions for the concentration-dependent col¬ 
lective diffusion coefficient and viscosity of permeable 
hard-sphere dispersions, which have been well tested 
against computer simulation results and experimental 
data on non-ionic PNIPAM microgels E3S3- These ex¬ 
pressions have been obtained using a hydrodynamic ra¬ 
dius model wherein the internal hydrodynamic particle 
structure is mapped on a hydrodynamic radius parameter 
for unchanged direct (non-hydrodynamic) interactions. 
Riest et al. have shown recently that this simple parti¬ 
cle model captures quantitatively the static and dynamic 
properties of non-ionic PNIPAM microgel dispersions in 
an organic solvent mi . This description is accurate 
up to large particle concentrations where the effect of sol¬ 
vent permeability is quite significant. Using a boundary 
layer analysis of the advective-diffusion transport in the 
CP layer US [M], we calculate the concentration profile 
and permeate flux at the membrane for different oper¬ 
ation conditions. Our focus is on particle permeability 
effects using concentration-dependent dispersion proper¬ 
ties. Various efficiency and energy cost indicators for the 
filtration process are introduced and numerically evalu¬ 
ated to check the efficiency and sustainability of the UF 
filtration process. 

The paper is structured as follows: in Sec. [TT] we 
describe the macroscopic cross-flow UF model and the 
similarity solution scheme for the CP boundary layer. 
We also discuss a transcendental equation that simplifies 
the problem for constant dispersion transport properties. 
Under such conditions, this equation leads to upper and 
lower bounds for the concentration profiles, which can 
be profitably used for a quick check on unwanted cake 
formation. The employed analytic expressions for the 
concentration-dependent collective diffusion coefficient, 
shear viscosity and osmotic compressibility of permeable 
Brownian particle dispersions are discussed in Sec. m 
Numerical results for the CP profile and the permeate 
flux are analyzed in Sec. |IV| regarding their dependence 
on the TMP, shear rate, feed concentration, particle size, 
and in particular the particle permeability. The signif¬ 
icance of the microgel permeability is assessed by com¬ 
parison with results for impermeable particles. We also 
evaluate the UF process in terms of process efficiency and 
cost indicators. Our conclusions are presented in Sec. [V] 

II. MACROSCOPIC UF MODEL 
A. Cross-flow transport 

We consider a homogeneous feed dispersion of non¬ 
ionic Brownian microgels particles of hard-core radius a 
dispersed in a Newtonian solvent at volume fraction <po- 
As illustrated in Fig. [l] the dispersion is longitudinally 
and steadily pumped through a hollow cylindrical fiber 
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permeate flux (solvent) 


FIG. 1. Sketch of inside-out cross-flow UF in a hollow cylin¬ 
drical fiber membrane. 


membrane of length L and inner radius R. Due to an ap¬ 
plied (transversal) transmembrane pressure, AP, a small 
fraction of the axially in-flowing solvent permeates the 
membrane to the outside of the fiber. We consider the 
membrane to be fully particle retaining. In this so-called 
inside-out cross-flow UF setup, the particle advection to¬ 
wards the membrane is balanced at steady-state by the 
diffusive back transport away from the membrane. 

A coarse-grained continuum mechanics description can 
be used for R a where the granular nature of the col¬ 
loidal Brownian particles, assumed to be monodisperse 
and neutrally buoyant, and the morphology of the mem¬ 
brane pores are not resolved. The momentum balance 
of the dispersion-averaged flow is then described by the 
effective Navier-Stokes equation, 

p + v • = - Vp + ?)Av 

+ Vr? • [Vv + (Vv) T ] . (1) 

We have used that the dispersion is incompressible 
on the macroscopic scale, implying that V • v = 0. 
In Eq. (JTJ), p is the constant dispersion mass density 
and r](<f>) is the effective dispersion shear viscosity de¬ 
pending on the particle volume fraction cj> = (47r/3) no 3 , 
with n denoting the mean particle number density. The 
dispersion-averaged velocity, v(r, t), and pressure, p( r, t), 
depend on the spatial position, r, and time t. Here, 
Vp = Vpf + Vn, wherein pp is the fluid-phase pres¬ 
sure contribution whose gradient adjusts itself such that 
the flow remains incompressible throughout the disper¬ 
sion, and n is the equilibrium osmotic pressure due to the 
Brownian particulate phase [22]. The third term on the 
right-hand side of Eq. |l]), invoking Vp and the dyadic 
tensor gradient Vv plus its transpose (Vv) T , originates 
from the concentration-dependent viscosity, i]((f>), in the 
dispersion viscous stress tensor. This term contributes in 
an inhomogeneous dispersion region such as the CP layer 
in UF. 

The mass balance (particle conservation) is described 


by the continuity equation, 


dn _ , 

di + v- J —o, 


( 2 ) 


where n(r,t) and J(r, i) are the coarse-grained local par¬ 
ticle number concentration and advective-diffusive parti¬ 
cle flux, respectively. The latter is the sum, 

J = J ad + J D =nv-DVn, (3) 


of advection and Brownian diffusion fluxes. Here, D(cf >) is 
the concentration-dependent long-time collective or gra¬ 
dient diffusion coefficient. For non-zero particle concen¬ 
tration and dominantly repulsive interactions, D(</>) is in 
general larger than the single-particle translational diffu¬ 
sion coefficient, Dq 1 with the latter measured at infinite 
dilution. We note that 

J« = -^)Vn = -Mvn, (4) 

is valid for a monodisperse dispersion of electroneutral 
particles at temperature T, with ks denoting Boltz¬ 
mann’s constant. Thus, the diffusion flux is driven by 
the gradient in the osmotic pressure n(</>), which in turn 
is proportional to the gradient in the particles chemical 
potential. The transport property K{<f>) is the long-time 
macroscopic sedimentation coefficient. It is defined as 
the ratio of the mean sedimentation velocity (measured 
relative to the center of volume rest frame) of a weakly 
settling homogeneous dispersion of volume fraction (f> to 
the settling velocity of an isolated particle. 

Substituting the particle flux into Eq. (pi) and using 
the local volume fraction rather than n( r, tp the macro¬ 
scopic advection-diffusion equation, 

|+v.V^V.(flV^), (5) 

is obtained. Analytic expressions for the collective diffu¬ 
sion coefficient, D(q i), and the zero-frequency dispersion 
viscosity, ?)(</>), of solvent-permeable Brownian particles 
are discussed in Sec. Mil 

Eqs. ([l])-([5]) apply to UF where Brownian parti¬ 
cle dispersions are considered under laminar flow condi¬ 
tions for which the single-particle shear Peclet number, 
Pe“ = rg, /P <0 . 01 , is so small that the system is only 
slightly perturbed away from thermal equilibrium (linear 
response regime). Non-Newtonian effects, such as normal 
stress differences and shear-thinning or thickening, are 
then absent or negligibly small. Here, rg> = a 2 /Dq is the 
Brownian diffusion time for diffusion across the particle 
radius and 77 = 1 /y is the shear flow advection time for 
the characteristic shear rate 7 . Assuming fully developed 
Poiseuille feed flow at the cylindrical fiber entrance, con¬ 
sistent with no-slip hydrodynamic boundary condition at 
the membrane surface, the characteristic shear rate is [23j 

_. _ 2 u m _ Ap L R 
~ 2r 1 (</> 0 )L' 


( 6 ) 
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where Ais the longitudinally applied pressure differ¬ 
ence across the fiber of length L , to be distinguished from 
the TMP, and u m is the inflow dispersion velocity at the 
fiber center (axis) equal to twice the cross-section aver¬ 
aged velocity (u). 

In typical microfiltration experiments Pe“ > 1 and 
Brownian diffusion is hence not important. The flow- 
driven particle transport at larger concentrations is dom¬ 
inated instead by non-equilibrium many-particle hy¬ 
drodynamic interaction (HI) effects, most notably by 
anisotropic collective hydrodynamic diffusion and cross¬ 
streamlines particle migration. These HI effects must be 
properly taken into account in the expressions for the 
coarse-grained particle current and dispersion-averaged 
stress tensor. The latter includes now normal stress dif¬ 
ferences, a shear-rate dependent non-equilibrium effec¬ 
tive viscosity, and an effective osmotic-type pressure con¬ 
tribution 181 31] , all depending on the local shear rate in 
addition to <j. i. Microfiltration is not treated in this work. 


B. Steady-state boundary layer 

Cross-flow UF is commonly performed at laminar flow 
conditions where the Reynolds number associated with 
cylindrical fiber flow fulfills that Re r = R{u m /v) < 
2000, with v = r)/p denoting the kinematic viscosity. 
This implies Re a = ( a/R)ReR -C 1 for the particle 
Reynolds number in terms of the bulk shear rate, guaran¬ 
teeing that particle inertia is negligible. Moreover, and 
quite importantly, the permeate velocity is small com¬ 
pared with the longitudinal cross-flow velocity in the bulk 
of the fiber, u m v%, where u„, is the maximal per¬ 
meate velocity measured for a clean membrane at given 
TMP using pure solvent as feed. Due to the strong ax¬ 
ial advection, and after a short transient [33 03, a sta¬ 
tionary thin CP layer develops on the membrane surface 
where diffusive-advective particle transport takes place. 
The thickness, S, of the diffuse boundary layer is much 
smaller than the inner radius, R 1 of the hollow fiber. This 
layer is characterized by marked axial concentration gra¬ 
dients. In the bulk region outside the boundary layer, 
the particles are transported essentially by axial advec¬ 
tion with preserved uniform concentration <j >o and unper¬ 
turbed flow profile. 

The locally planar CP boundary layer is conveniently 
analyzed using the local Cartesian coordinates shown in 
Fig-! where the x-axis with the unit vector e x extends 
longitudinally along the membrane surface and the y-axis 
with the unit vector e y stretches out radially towards the 
fiber axis. The steady-state flow velocity in this local 
frame, 


Using d/dx ~ 1/L and d/dy ~ 1/5 with 6 <C L in Eq. 
([5]), the steady-state advection-diffusion equation reduces 
in the CP layer to 

U[X ' V) f I +V{X ' V) fy = Ty{ DW fy)' (9) 

where the axial diffusion flux contribution has been ne¬ 
glected in comparison to the transversal one. A similar 
analysis largely reduces the effective Navier-Stokes Eq. 
(|T|) in the CP layer to the remaining expression, 



describing constant axial dispersion shear stress across 
the CP layer. The stress is independent of the axial posi¬ 
tion x and is obtained from the asymptotic matching to 
the bulk region as [23,333 


v(4>) 


du 

dy 




y —± 00 




(ii) 


Note here that the time scale for which steady-state 
viscous flow has developed is much shorter than the 
time scale for the development of the steady-state CP 
layer. This is due to the very large Schmidt number, 
Sc = v/D 0l in colloidal dispersions, given by the ratio of 
fluid vorticity and particle thermal diffusion coefficients 

@12051. 

The steady-state and positive valued permeate velocity 
through the uniformly porous membrane with its inner 
wall at x = 0, 


v w {x) = -v(x), 

is described in UF using Darcy’s law by 

n (<t> w {x))' 


>{x) = 


1 - 


A P 


( 12 ) 


(13) 


where 


v° w = L p AP 


(14) 



v(r) = u(x, y) e x + v(x, y) e y , 


fulfills the incompressibility condition 


du dv 
dx dy 


(7) 

( 8 ) 


FIG. 2. Locally flat CP layer next to the surface, y = 0, of 
the hollow cylindrical fiber membrane of length L in steady- 
state cross-flow UF. The osmotic-pressure-regulated permeate 
transverse flow profile, v{x, 0) = — v w (x), and the axial bulk 
flow profile, u(x,y), are sketched. 
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is the permeate velocity at zero particle concentration, 
for a given TMP A P > 0. Here, L p = l/(r]oR m ) is the 
hydraulic permeability of the clean membrane, and R m 
is the associated hydraulic resistance. According to Eq. 
(13), the permeate velocity is determined by the differ¬ 
ence between TMP and the osmotic pressure, n(^ UJ ), at 
the inner membrane wall. Here, 


4> w (x) = <j)(x, 0 ) 


(15) 


is the particle volume fraction at the inner membrane sur¬ 
face which increases with increasing distance x from the 
fiber entrance. In obtaining Eq. (13), we have assumed 


(i) that the membrane is impermeable to the particles, 

(ii) zero colloid concentration in the permeate, (iii) sta¬ 
tionary, fully developed permeate flow inside the mem¬ 
brane, and (iv) the absence of fouling effects which ef¬ 
fectively enlarge the hydraulic resistance above the value 
R m of a clean membrane. We also neglect the axial pres¬ 
sure drop which in typical UF applications is far smaller 
than the TMP. 

The remaining boundary conditions used in solving 
Eqs. ([ 8 ])-( 10) are as follows: we use the no-slip boundary 


condition for the axial velocity component at the mem¬ 
brane surface, 


u(x, 0 ) = 0 , (16) 

and assume that the particle concentration right at the 
fiber entrance and distant from the membrane surface 
matches the bulk (feed) concentration, 


C. Similarity solution 


The similarity solution scheme used in this work fol¬ 
lows that by Shen and Probstein for gel polarization [23]. 
We account here in addition for the osmotic pressure con¬ 
tribution to UF and give additional insight into the ob¬ 
tained solution. 

We first estimate the thickness, 8 = 8(x), of the CP 
layer. To this end, we use that the transverse advection 
term, vd <j>/dy, in Eq. ([9]) is negligible at the outer dif¬ 
fuse edge of this layer in comparison to the axial one, 
since u m and since <f> has decayed basically to ipQ. 

Right at the membrane surface, however, the axial ad¬ 
vection term is zero instead due to the no-slip boundary 
condition. Using that d/dy ~ l/<5, the estimate 


8 = 



(19) 


is obtained, with Db = D((f>o) denoting the bulk collective 
diffusion coefficient. Eq. (19) describes an axial growth 
of the CP layer thickness proportional to a; 1 / 3 . It follows 
that 8(L)/R ~ (L/R) 1 / 3 / (Pe^) 1 ^ 3 . The axial advection 
Peclet number Pe^ = r^/rad, with = R 2 /D 0 and 
Tad = R/u m , is in UF orders of magnitude larger than 
one, consistent with the thin boundary layer assumption 
we have started from. 


Eq. (19) implies the existence of a similarity solution, 


</>(A), for the CP profile which for each axial position 
x is a function of the dimensionless similarity variable 


</>(0,y) = 4>(x,y oo) = </>o. (17) 



( 20 ) 


Finally, the assumed impermeability of the membrane 
regarding the colloidal particles is considered using the 
zero normal particle flux condition, 


only with A oc y/x 1 ^ 3 . The dispersion incompressibility 
condition in Eq. ([ 8 ]) and the shear stress Eq. (11) are 
expressed in terms of A as 




■V w (x) <j) w {x) = 0, 


(18) 


y— o 


at the membrane surface. In using this reflecting bound¬ 
ary condition, particle adsorption by the membrane is 
excluded which adds to fouling. 

For given dispersion properties D(</>), and n(</>), 

the non-linear coupled partial differential Eqs. (fs])-( fTTTI) i n 
conjunction with the boundary conditions in Eqs. ( |ll| )- 
(18) form a two-dimensional boundary value problem. It 


can be solved only numerically for the searched for CP 
concentration and velocity profiles <j>(x,y) and v(x, y), 
respectively. The only exception is the limiting case, 
= 0 , of an impermeable fiber wall where <f> = </> 0 
holds uniformly throughout the fiber and the Poiseuille 
velocity profile extends right up to the membrane wall. 
The solution of the boundary value problem for > 0 
is largely simplified using a similarity analysis described 
in the following. 


dv 8 ^du 

dX 3x dX 

1 du 1 
j8 dX 1 

respectively, where f)( 0 ) = ri{4>) / f]{4>o) is the reduced vis¬ 
cosity. 

Dimensionless velocities are introduced accordingly by 


( 21 ) 

( 22 ) 


U(A) 

U{x,y) 


3x 

~~^ v 

u 

j8 1 


(23) 

(24) 


where the axial velocity u has been scaled with the fluid 
velocity, 7 <5, at the outer edge of the CP layer, and the 
transverse velocity v by the velocity D^/S = S 2 j/(3x) 
characteristic of the advection-diffusion balance. The 
negative sign in the definition of V has been included 
to make it positive valued. Different from V, the re¬ 
duced axial velocity U is not a function of A only. This 
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follows from considering the special case, 77 = 1 , of con¬ 
stant viscosity where u = 7 y is the solution of Eq. (22) 
satisfying the no-slip condition at the membrane sur¬ 
face. The corresponding transverse velocity is obtained 
from the incompressibility condition and Eq. (12) as 
v(x,y) = —v w (x), being constant across the CP layer 
for a given x. 

The non-dimensional velocity related to U defined as 


must be calculated for each x value separately. We have 
performed these calculations with MATLAB’s routine 
bvp4c 071, using analytic expressions for collective diffu¬ 
sion coefficient, effective dispersion viscosity and osmotic 
pressure presented in the following section. 


The general behavior of the concentration profile fol¬ 
lows from the integrated Eq. (28), 


G(A) = U — X, (25) 


is, however, an explicit function of A. Substitution of Eq. 
(l25|) into Eq. (p2l) leads to 


dG _ 1 

dX ) 


(26) 


which is the non-dimensionalized shear stress equation. 
The dimensionless incompressibility condition follows 
from Eqs. (21) and (25) as 


dV 

~dX 


= G(X) - X 


dG 

~dX' 


(27) 


= _^exp{_ f X dX'm±^l}, 
dX D(<j>(X)) l J 0 D(</>( A')) J 

(35) 

where the zero flux condition at the membrane wall has 
been used. In accord with physical expectation, it de¬ 
scribes the strictly monotonic compressed exponential 
decay of <j)(r ]) from the wall value <f> w at A = 0 towards 
the large-A bulk value (f> 0 where dtp/dX —» 0. 


D. Solution for constant viscosity and diffusivity 


Using Eqs. 120J) ([24]), ([25]) and ([27]), the advection- 
diffusion Eq. Q is expressed in terms of A, V and G 
as 


d_ 

dX 




= 


dX l 


A [A + G(A)] + V(A) j , (28) 


where D(<p) = D^/D^. It is noticed here that with q i 
and V, also G is a function of A. 

In terms of A, V and G, the boundary conditions in 
Eqs. ([l2]), dl 6 ])-([l 8 ]) are expressed as 


G( X = 0) = 0 
V(A = 0) = V W 
4>{x -» 00 ) = 4 > 0 


dcf) 


D{4>w)~7t (A — 0) + V w <f> w — 0 . 
dX 


(29) 

(30) 

(31) 

(32) 


Here, V w = ( 3x/6 2/ y)v w is the dimensionless permeate 
velocity, with v w given by Eq. (13). Using the boundary 


conditions at A = 0, Eqs. (26) and (27) are integrated, 


G(A) = / 

Jo 


dX' 


0 r/W’O')) 


V(A) - V w = 2 [ dX 1 G(X') - XG(X ), 
Jo 


(33) 


(34) 


with G and V expressed now as functionals of the con¬ 
centration profile. 

Eqs. (26)-(28) form a set of coupled nonlinear ordi¬ 


nary differential equations which in combination with the 
boundary conditions in Eqs. (29)-(32), can be solved nu¬ 


merically for the searched for profile <fi(X) and the reduced 
velocities G(A) and U(A). Since V w in Eq. (13) varies 


with the axial position along membrane surface, due to 
the likewise varying osmotic pressure, the set {tp, G, V} 


An additional simplification follows for a 
concentration-independent reduced viscosity fj in 
the CP layer not necessarily equal to the bulk value of 
one. In this case, the solutions of Eqs. ([26]) and 27) are 


G = A(l /77 — 1 ) and V = V w , respectively, for all A > 0. 
Assuming in addition a cons tant reduced collective 
diffusion coefficient D, Eq. (35) can be integrated using 
Eq. (31), resulting in [21] 


<t>o 


1 - (V w /D) /“ dAexp{—A 3 /(3t)D) - A V w /D} ’ 

(36) 

where V w = V w {(f > u ,). This is a transcendental equation 
for the membrane wall concentration, <fi u ,(x), at a given 
axial distance x from the feed inlet. It can be solved us¬ 
ing an iteration method. For given osmotic pressure ex¬ 
pression, the permeate velocity v w {x ) follows then from 


substituting the solution <j) w (x) into Eq. (13). In Sec. IV 


Eq. (36) is profitably used to generate a good approx¬ 
imation of the profile <f> w (x) for non-constant collective 
diffusion and viscosity coefficients. 


III. PERMEABLE PARTICLES DISPERSIONS 


We model non-ionic microgel dispersions generically 
as dispersions of uniformly solvent-permeable colloidal 
hard spheres [571 - HD] . As shown in a recent theoretical- 
experimental work by Riest et al. m, this simplify¬ 
ing model provides an excellent description of static and 
dynamic light scattering data, and hence the pair dis¬ 
tribution function and collective diffusion properties, of 
non-ionic PNiPAM microgels m ■ 
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A. Osmotic pressure 


The Brownian particles accumulated in the CP layer 
give rise, for Pe“ <C 1, to an equilibrium osmotic pressure 
profile II(</>). According to Eq. Q, the osmotic pressure 
gradient drives a thermal diffusion current pointing away 
from the membrane into the bulk. On the other hand, 
the osmotic pressure at the membrane wall, re¬ 

duces the permeate flux as quantified by Eq. (13). Since 


the microgels are modeled as permeable hard spheres, 
which are characterized statically by the hard-core ra¬ 
dius a with associated volume fraction <j>, the osmotic 
pressure in the fluid phase is described to excellent ac¬ 
curacy by the Carnahan-Starling (CS) equation of state 


m = 


i+ cj)+(f> 2 -1 

mF 


(37) 


where Z((j>) = 11(0)/ (n(</>)fcsT) is defined as the ratio of 
osmotic and ideal gas pressures. Eq. (37) applies not only 


to the fluid-phase branch extending up to the freezing 
concentration, <j>f ~ 0.494, but decently well also to part 
of the non-crystalline metastable phase branch from cj>f 
up to the melting volume fraction cj) m ss 0.545. The CS 
Z{4>) does not reproduce the osmotic pressure divergence, 
II ~ 1/(0 — 0 rcp ), of the metastable non-crystalline phase 
at the random close packing concentration <f>- 


rep 


0.64 

This divergence is of no concern in this work 
where stagnant cake layers formed by particle jamming 
are not considered. 


B. Hydrodynamic particle model 

The nresh-size averaged flow inside a microgel particle 
is described in our model by the Brinkman-Debye-Bueche 
(BDB) equation (49j [5Q3, wherein the amount of fluid 
permeability is characterized by the non-dimensional pa¬ 
rameter 


X = Ka, (38) 

equal to the inverse of the hydrodynamic penetration 
depth, ft -1 , in units of the hard-core radius a. The pen¬ 
etration depth is equal to the square root of the Darcy 
permeability of a microgel particle, and it roughly equals 
its mean mesh size. Values of \ can vary in principle 
between x = oo, corresponding to an impermeable hard 
sphere with no-slip surface boundary condition, down to 
low values characteristic of strongly permeable particles. 
For the continuum picture underlying the BDB equation 
to be valid, the mean pore size should be no larger than 
one tenth of the particle radius (x > 10). Typical values 
in experimentally studied dispersions are \ ~ 20 — 40 or 
larger [40) [511. 

Transport properties of permeable particle dispersions 
such as D and depend on \ in addition to <j>. At infinite 


dilution, D(\, 4>) reduces to the single-particle transla¬ 
tional diffusion coefficient of a permeable sphere, 


D 0 (X) 


k B T 

6nr] 0 ah(x) ’ 


(39) 


with the hydrodynamic particle radius, ah, given by E3 


ah(x) = 2y 2 [ x - tanh(y)] 

a 2 X 3 + 3 [y - tanh(x)] ' 


(40) 


The hydrodynamic radius becomes equal to a for x -A 
oo. It decreases with increasing permeability (decreasing 
y) since the viscous drag is reduced when more fluid is 
allowed to penetrate the particle. For realistic param¬ 
eter values, x ^ 10, the particle-fluid interface can be 
considered as locally flat. The ratio ah/a is then well 
approximated by the flat interface value 

= *-(■ (41) 


a hj(x) 


with the error introduced in ah being of order 0( 1/x 2 ) 
small. The flat-interface radius, ahj, and the associated 
reduced slip length, L^j = 1 — ahj/a = 1/x, are inde¬ 
pendent of the single-particle transport property used for 
their definition. In other words, the same value for Lhj 
is obtained regardless whether the translational or rota¬ 
tional single-particle diffusion coefficients or the intrinsic 
viscosity are used mm- 

As explained in Refs. ESI S3, diffusion and rheologi¬ 
cal properties of a large variety of dispersions of spher¬ 
ical colloidal particles of quite different internal hydro- 
dynamic structure and surface boundary conditions, in¬ 
cluding core-shell particles and particles with Navier slip 
boundary conditions, can be characterized hydrodynam- 
ically in terms of the single hydrodynamic parameter ah- 
This key observation allows for the usage of the simpli¬ 
fying hydrodynamic radius model wherein the particles 
are described hydrodynamically as no-slip spheres, char¬ 
acterized by ah derived from a single-particle transport 
property, for unchanged direct particle interactions. Re¬ 
garding the considered non-ionic microgels, this means an 
unchanged direct interaction radius a. The error intro¬ 
duced in using the hydrodynamic radius model for trans¬ 
port coefficient calculations, such as those for D a nd y 

is of O ({L* h y) 2 ^ small [ID]. Since Eqs. (401 and (41) 


describe a one-to-one mapping between x and the uni¬ 
versal parameter ah, the at first sight specific model of 
uniformly permeable spheres used in this work is in fact 
quite general. 


C. Collective diffusion 

For the considered macroscopic length and time scales, 
collective thermal diffusion or gradient diffusion refer to 
the relaxation of thermally induced local concentration 











gradients by the cooperative diffusive motion of Brown¬ 
ian particles which interact both directly and hydrody- 
namically. The long-time collective diffusion coefficient 
appearing in Fick’s law of macroscopic gradient diffu¬ 
sion of a non-ionic dispersion of monodisperse permeable 
spheres is given, according to Eq. @,by 

D{xA) = D 0 {x) Ki ^, (42) 

namely by the single-particle diffusion coefficient Dq(x) 
times the ratio of the long-time sedimentation coefficient, 
K(x, 0) and the thermodynamic coefficient S(<f>). The 
latter is equal to the zero-scattering wavenumber limit 
of the static structure factor determined in scattering 
experiments [52| and to the reduced isothermal osmotic 
compressibility, 


i ( amm 

V <90 ) ^ gT 


(43) 


where /i s denotes the chemical potential of the dispersing 
fluid. 

An issue not addressed in the UF literature is that the 
long-time collective diffusion coefficient D , and the as¬ 
sociated long-time sedimentation coefficient K , in Eqs. 
© and ( |42[ ) are smaller than their respective short-time 
counterparts. This is due to the, at long times, slow¬ 
ing influence of the micro-structural environment, which 
is perturbed from the equilibrium isotropic state by the 
sedimenting or collectively diffusing particles [53j. Dif¬ 
ferent from self-diffusion and viscosity, where the respec¬ 
tive short-time and long-time transport coefficients are 
strongly different in concentrated dispersions, the short- 
and long-time collective diffusion and sedimentation coef¬ 
ficients differ by less than 7% only, even in concentrated 
dispersion of no-slip hard spheres where the hydrody¬ 
namic interactions are particularly strong 1541 . The rel¬ 
ative difference in the latter two quantities is of purely 
hydrodynamic origin and caused by the non-pairwise ad¬ 
ditive near-distance part of the hydrodynamic particle 
interactions [ 55] , We emphasize that the present discus¬ 
sion applies to Brownian particles for low sedimentation 
Peclet number characteristic of UF. In the large Peclet 
number regime, micro-structural changes and their effect 
on the concentration dependence of the sedimentation 
velocity become highly significant [56] . 

In this work, we use the short-time form as a substitute 


for the long-time K(x, ( t ) ) in Eq. (421 since first, the cal¬ 


culation of the short-time sedimentation coefficient is far 
simpler and second, precise simulation data are available 
for the short-time coefficient in its dependence both on 
0 and % [S3 SSI HD]- The wall concentrations in our UF 
results are such that (f> w < 0.4 so that the long-time coef¬ 
ficients D and K are only slightly overestimated by this 
substitution. The inset of Fig. 3 shows simulation re¬ 
sults for the concentration dependence of the short-time 
K(x, (j>), for values of the inverse reduced permeability 



<l> 


FIG. 3. Inset: concentration dependence of the short-time 
sedimentation coefficient, K(x, 0), of permeable hard-sphere 
dispersions. Main figure: corresponding reduced short-time 
collective diffusion coefficient, D(x,4’)/Do{x), obtained from 
iv(x, 0) by division with the CS S(0). O and o are simula¬ 
tion results for impermeable spheres (x = oo) by Ladd [57] 
and Abade et al. 22- 1=1 j A, and V are simulation results 
by Abade et al. [37] for x = 100, 50, and 20, respectively. 
Colored lines are polynomial fits to the respective simulation 
data. Dotted and dashed lines: Eq. ( |46| ) for a = 6.55 and 
the cell model prediction by Happel 1341 . respectively, both 
for impermeable spheres. The cell model predictions for K 
and D are non-analytic at 0 = 0. Black solid lines in main 
figure and inset: second-order virial expan sion results for im¬ 
permeable spheres given in Eqs. (451 and (44 1 , respectively. 


given by x = 20, 50, 100, and oo. These values corre¬ 
spond, according to Eq. (41), to reduced hydrodynamic 


radius ah/a = 0.95, 0.98, 0.99, and 1. 

For impermeable no-slip hard spheres, the numerically 
precise second-order virial (concentration) expansion re¬ 
sult, 


A'(0) = 1 - 6.5460 + 21.9180 2 + O (0 3 ), 


(44) 


for the short-time sedimentation coefficient has been ob¬ 
tained by Cichocki et al. [58]. As seen from the inset 
of Fig. [3j this expression describes the simulation data 
[S3 [57] well for 0 < 0.1 only (solid line). The division 
of K (0) in Eq. (44) by the second-order virial expansion 
result of the hard-sphere osmotic compressibility, 
leads to the second-order virial expansion expression [59] 


D{<j>)/D 0 = 1 + 1.4540 - O.450 2 + 0(0 3 ), 


(45) 


for the short-time collective diffusion coefficient. This 
quadratic polynomial describes the weakly monotonic in¬ 
crease of D(0) with increasing 0 astonishingly well even 
up to 0 « 0.5 where the non-sheared dispersion starts to 
solidify. See here the solid black line in the main part of 
Fig. § This finding points to a mutual cancellation of 
higher-order terms in the virial expansions of K and S. 
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The simulation data for D shown in the main part of the 
figure have been obtained from dividing the simulation 
data for K by the CS 5(0) obtained from Eqs. ([37]) and 

©• 

Previous UF works dealing with impermeable colloidal 
particles have commonly used either simply a constant 
D = Do , or the semi-empirical expression GEEim 


KM) = (1 - <t>T 


(46) 


with the exponent a = 6.55 selected to match Batchelor’s 
[53] first-order virial coefficient result (Eq. (44) and the 
dotted line), or the cell model result for K( 0) by Hap- 
pel |34] (dashed line in inset) in conjunction with the CS 
S((j)) as divisor [2EHS3- While Eq. ( |46| ) describes the 
simulation data for impermeable spheres depicted in the 
inset overall decently well, the decay of AT (0) at larger 0 
is overestimated to such an extent that, after division by 
5(0), the resulting D is decreasing for larger 0, in con¬ 
flict with simulation data and experimental findings for 
colloidal hard spheres |59j . Notice here that K(cf >) and 
S((j>) are both monotonically decreasing with increasing 
0 in a concave-shaped form. Their ratio is thus very 
sensitive to small changes in both quantities. Using the 
exponential form for K in Eq. (46), Bhattacharjee et al. 


EQ] combined an approximate long-distance treatment of 
the hydrodynamic mobilities of two no-slip spheres with 
the Percus-Yevick solution for the pair distribution func¬ 
tion leading to a concentration-dependent exponent cn( 0 ). 
While their modified approach predicts the monotonic 
increase of DM) with increasing 0 , the increase is sig¬ 
nificantly overestimated in comparison to the simulation 
data. 

The cell model predictions for I\ and D (dashed lines) 
are ruled out since they violate the exact analytic low -0 
forms of these quantities noted in Eqs. (44) and (45), 


respectively. In particular, an infinite negative slope of 
DM) at 0 —* 0 is predicted by this model. In summariz¬ 
ing the present discussion on collective diffusion of imper¬ 
meable particles, our conclusion is that Eq. (451 should 


be used as a good description of the DM) of impermeable 
colloidal hard spheres. 

A thorough analysis of the sedimentation coefficient 
of colloidal dispersions of hard spheres with internal hy¬ 
drodynamic structure (non-zero permeability) was given 
in [40j . This article provides an approximate analytic 
expression for A'(y, 0), where 7 = ah/a, describing the 
simulation data quite well for all (j> < 0.5 and 7 > 0.8 
using Eq. (41). We only quote here the resulting low-0 
expression 


DjxM) 

A)(x) 


1+ 1.454 + 


8.592 

X 


0+0(0 2 ), (47) 


valid for y > 10 , with 0( 1 /y 2 ) corrections disregarded. 
It quantifies the initially linear increase of K(xM) at 
low 0 displayed in Fig. [3] with increasing slope with 
increasing permeability. 


The remaining smaller differences between the analytic 
expression for K in |40| and the simulation data are en¬ 
hanced at larger 0 after the division by 5(0). In our UF 
calculations for permeable particles, we therefore use a 
direct fit of the simulation data for D(x, 0) EH in the 
form of a fifth order polynomial (colored lines). 

The simulation data for K and D for a given concen¬ 
tration increase with increasing permeability (decreasing 
y). The respective increase is more pronounced for larger 
X since the main effect of permeability is to lower the hy¬ 
drodynamic interactions between the particles. Perme¬ 
able hard spheres sediment and diffuse collectively faster 
than impermeable ones for two related reasons: first, the 
hydrodynamic interactions are weaker and second, Dq 
is larger due to ah < a. The sedimentation coefficient 
decreases with increasing volume fraction owing to an 
increasing fluid backflow and enhanced hydrodynamic 
interactions El ISO]. For a given y, D increases with 
increasing 0 which reflects that the diminution in K is 
overcompensated by the diminution in the osmotic com¬ 
pressibility. The monotonic increase of D with increas¬ 
ing 0 is characteristic of colloidal hard spheres but not of 
particles with longer-ranged repulsive interactions or par¬ 
ticles having shorter-ranged attractions. For example, in 
lower-salinity dispersions of charge-stabilized colloids and 
protein solutions, D behaves non-monotonically, passing 
through a pronounced maximum with increasing 0 m- 


D. Effective dispersion viscosity 


The effective dispersion viscosity, rj , in a steadily 
sheared dispersion is the ratio between the dispersion 
shear stress and the rate of strain [35], It depends on 
the concentration and internal hydrodynamic structure 
of the Brownian particles and their direct and His. Since 
Pe“ < 1 in UF, there are no shear-thinning effects. The 
zero-frequency low-shear-rate limiting effective viscosity 
of a permeable particles dispersion is the sum 


viXi 0) = VooixM) 


1 A?/(y, 0 ) ~ 

Voo(xM)\ ’ 


(48) 


of the high-frequency viscosity contribution, 77 ^ (y, 0 ), 
which is of purely hydrodynamic origin, and the shear 
relaxation viscosity part, Arj(xM)i originating from the 
relaxation of the shear-perturbed dynamic particle cages 
formed around each particle. The viscosity part, A 77 , is 
influenced both by direct and His, with the consequence 
that, for concentrated dispersions, the long-time viscosity 
77 is substantially larger than the short-time viscosity 77 ^. 
This distinguishes the viscosity from collective diffusion 
where the difference between the associated short-time 
and long-time collective diffusion coefficients stays small. 


As discussed in Refs. 00] and 


. the term in brackets 


in Eq. (48) can be expected to be only weakly influenced 
by His. This suggests the no-HI factorization approxi- 
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mation 


&v{xA) _ ^v{4>) 


Voo{xA) Vo 


(49) 


D —HI 


where A(?^)| no _Hi is the shear relaxation part without HI, 
with tJqq | no —hi being equal to the fluid viscosity 770 - In 
this approximation, the hydrodynamic particle structure 
quantified by y and the His are assumed to affect 77 only 
through the factored-out high-frequency viscosity in Eq. 
(48 1 . This simplifies the calculation of A 77 considerably. 


An analytic estimate for the A 77 of colloidal hard spheres 
without His is given by 


At?(0) 


Vo 


D-m 


^<)> 2 5(2a + , 4>), 

5 


(50) 


where g(2a + ,cf>) = — l)/4</ is the contact value 

of the equilibrium radial-distribution function, and Z((f>) 
is given for <j> < 0.5 by Eq. (37). This estimate 


combines the exact quadratic-order low concentration 
limit of A(? 7 )| no _Hi with its divergence at the random 
closed packing volume fraction, 0 rcp , constituting the up¬ 
per concentration limit of the non-crystalline metastable 
branch of the hard-sphere phase diagram. This diver¬ 
gence is triggered by the divergence of g(2a + ,cj)) at <^> rcp 


The existing high-precision simulation data |38] for 
the high-frequency viscosity of permeable colloidal hard 
spheres in the fluid phase region <f> < 0.5 are well de¬ 
scribed by the the generalized Saito formula ESI HU, 


Vooix,<t>) = l + [v\(x)</>- 


i + S(xA) 


1 ^ !</>M(x) 


1 + s(x,</>) 


, (51) 


with the Saito function S{x,4>) = '/>(& h(x ) — 2/5) [ 77 ](%). 
Both the intrinsic viscosity, given by [ 77 ] (y) ss 
(5/2) (1 — 3/y) for x > 10, and the Huggins coefficient, 
kh(x)> depend on the particle permeability parameter y 
only. Analytic fitting expressions of these two quantities 
are given in Ref. TO . 

In this work, we use the factorization approximation 
in Eqs. (48)-(51) as the analytic input for ?)(y, 0) = 
vix> 4^/viXi </>o7"ui our UF calculations. 

The concentration dependence of 77 (y, </) is plotted in 
the main part of Fig. [3] together with the associated high- 
frequency viscosity, v 00 (x.^). shown in the inset. While 
both viscosities increase with increasing concentration, 77 
is substantially larger than 77 ^ for large concentrations, 
due to the at large 4> dominating shear relaxation contri¬ 
bution A 77 . Notice the different scales of the ordinate in 
the inset and the main figure part. 

We consider first the high-frequency viscosity pre¬ 
sented in the inset. The simulation data by Abade et 
al. (open symbols) for permeable hard spheres {35] show 
that 7700 decreases with increasing particle permeability 
(decreasing y) since viscous dissipation is reduced. The 
simulation data are quantitatively described by the gen¬ 
eralized Saito formula in Eq. (51) for all considered per¬ 


meabilities (solid lines in the inset). For comparison, also 



<l> 


FIG. 4. Concentration dependence of the steady-shear zero- 
frequency dispersion viscosity ? 7 (x, 0 ) (main figure), and the 
associated high-frequency viscosity goo{\,(j>) (inset), for per¬ 
meable colloidal hard spheres. •, ■ and O are experimental 
data for 77 by Segre et al. [52], and Weiss et al. EE and 
simulation results by Foss and Brady [64], respectively, all for 
X = 00. Solid black, red, violet, and blue lin e s: fa ctorization 
approximation predictions for 77 in Eqs. (48 1 -( 511 at x = 00 , 
100, 50, and 20, respectively. Black dotted line: Krieger- 
Dougherty formula for x — 00 EH) using <)> max = <j>c P ■ In 
the inset: O and o are simulation data for 7700 by Foss and 
Brady [64] and Abade et al. ESI. for impermeable spheres 
(y = 00 ). □, A, and V are simulations results by Abade et 
al. EH] f° r X = 100, 50, and 20, respectively. Correspondingly 
colored solid and dashed lines in the inset are the results by 
the generalized Saito formula in Eq. (511, and the cell model 
expression by Ohshima [55], respectively. 


recent cell model results by Ohshima [B5] for the 7700 of 
permeable spheres are shown (dashed lines). For imper¬ 
meable particles, y = 00 , the cell model prediction repro¬ 
duces an earlier result by Ruiz-Reina et al. [ 66 ] . As it can 
be noticed from the inset, the cell model method strongly 
overestimates the high-frequency viscosity at larger <f> for 
all considered permeabilities. 


Consider next the concentration dependence of the 
steady-shear viscosity, 77 , displayed in the main part of 
Fig. 0 The experimental data (filled black symbols, 
taken from [621163] ) and in particular the simulation data 
(open diamonds, from |64j ) for impermeable hard spheres 
are overall well represented by the factorization approx¬ 


imation in Eqs. (48 )-rt50h combined with the generalized 
Saito formula in Eq. ( |51[ ) for 77 ^( 00 , </>) (solid black line). 
A slight improvement over this factorization approxima¬ 
tion result for 77 has been discussed by Riest et al. J40] , 
This improvement, however, does not noticeably change 
our UF predictions for the CP layer and permeate flux 
profiles. 

Analogous to the corresponding behavior of 77 ^ (y, <f >), 
the steady-shear viscosity 77 is lowered with increasing 
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permeability due to the weakened His. This trend is 
described by the factorization approximation results for 
the three considered finite values of \ (differently colored 
solid lines). To our knowledge, no simulation or exper¬ 
imental zero-frequency viscosity data of permeable par¬ 
ticles are available to date for a direct comparison with 
our theoretical predictions for r/(y, </>). 

In previous UF works dealing with hard-sphere-like 
dispersions |2ftll32j . either a concentration-independent 
value for 77 or the phenomenological expressions by Fil¬ 
ers or Krieger-Dougherty have been used. The Krieger- 
Dougherty expression for the steady-shear viscosity reads 

m 


JL 

Vo 



(52) 



where the particle permeability enters only via the intrin¬ 
sic viscosity in the exponent. Here, </> max is the maximal 
volume fraction at which the expression diverges. In pre¬ 
vious UF works ED 131, the closed packing fee crystal 
value 4> cp = \/27t/ 6 ~ 0.74 was used for </ max . According 
to the dotted line in Fig. [4] valid for \ = 00, the Krieger- 
Dougherty relation for </ max = <fi cp strongly underesti¬ 
mates the steady-shear viscosity at large concentrations. 
This underestimation is less severe if, as suggested by de 
Kruif et al. m for low shear Peclet numbers, </> rcp ss 0.64 
is used instead of </> cp . The Eilers expression describes a 
viscosity curve (not shown in the figure) similar to the 
one representing the Krieger-Dougherty equation. 


IV. FILTRATION MODEL RESULTS 


FIG. 5. Concentration profile, cj) w {x ), of impermeable col¬ 
loidal hard spheres at the membrane surface, as function of 
the axial distance, x, from the inlet, in units of the fiber 
length L. Black solid line: similarity solution scheme re¬ 
sult for conc entration-d ependent and 77 ( 0 ) according 

to Eqs. (451 and Eqs. (481-(511, respectively. Red dashed 
and blue dotted lines: similarity solution results for constant 
D((fr 0 ) ~ Do and r]{(j> 0 ) as 770 , and for variable D(<)>) and con¬ 
stant 0 ), respectively. In the inset, the similarity solution 
result for concentration-variable D and 77 (black solid line) is 
compared with the solutions of the transcendental Eq. < |36| ) 
for constant D{(j> w ) and 77 (f)) 0 ) (red dashed-dotted line), and 
constant D{ij> 0 ) and ?7(</>u,) (blue dashed-dotted line). The 
black dashed-dotted line is the arithmetic average of the lat¬ 
ter two solutions. System parameters: A P = 5 kPa, 7 = 65 
s _1 , a = 10 nm, (j >0 = 1.0 x 10 
mPa - 1 s -1 . 


and L p = 6.7 x 1CT 10 


We present in this section our boundary layer model 
results for cross-flow UF of permeable (more generally: 
hydrodynamically structured) Brownian particles disper¬ 
sions. The selected system parameters are: membrane 
length L = 0.5 m, inner radius R = 0.5 mm, hydraulic 
membrane permeability = 0.5—6.7x 10 -10 mPa - 1 s -1 , 
and feed particle volume fraction </>0 = 1.0 — 1.5 x 10 -3 
which is small enough that the Einstein expression 77 ss 
(5/2) (ah/a) 3 4>o for permeable particles, with ah/a ac¬ 
cording to Eq. (40), applies for the feed viscosity. The 
calculations are performed for the room temperature 
value T = 293.15 K. 

The explored parameters in our calculations are the 
applied TMP A P = 0.3 — 5.0 kPa, particle excluded vol¬ 
ume radius a = 10 — 30 nm, and characteristic shear rate 
7 = 30 — 130 s -1 , in addition to the inverse permeability 
parameter \ = 20, 30, 50, 100, and 00 . These are realistic 
operation conditions for the UF of permeable Brownian 
particles, which we have selected such that Pe“ <C 1, 
with Brownian diffusion dominating shear-induced hy¬ 
drodynamic diffusion [T2]. Moreover, since R <C L, the 
boundary layer condition u m and, with Rea ~ 1 , 

also the condition of laminar pipe flow are fulfilled. The 
pipe radius R and the feed velocity u m enter only implic¬ 
itly into the similarity solution scheme for the CP layer 


through Eq. ([ 6 ]) for the characteristic shear rate. The 
particle radius a affects the CP layer and permeate flux 
directly by the factor D 0 (%) in Eq. (42). We further 
note that the selected operating conditions are such that 
membrane fouling caused by particle jamming is avoided. 
Fouling is not considered in the present UF study which 
focuses on particle permeability effects on the CP layer 
and permeate flux, based on accurate expressions for the 
concentration and permeability dependent viscosity and 
collective diffusion coefficient. 


A. Concentration-dependent transport properties 

Before discussing the influence of the particles perme¬ 
ability on cross-flow filtration, we analyze the effects of 
the concentration dependence of the transport proper¬ 
ties and 77(f))). For later comparison, we consider in 
this subsection the limiting case, \ = 00, of impermeable 
hard spheres. 

Figure [5] depicts the CP layer profile at the membrane 
surface, <p w (x ), in its dependence on the reduced axial 
distance, x/L, from the fiber inlet (Fig. [ 2 ]). The main fig- 
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lire part includes results of the similarity solution scheme 
discussed in Subsec. Eel The black solid line is the re¬ 


sult for (j) w {x) based on Eq. (45) for D(4>)/D 0 , and Eqs. 
(48)-(|5T|) for r]((j))/ri 0 . It shows the typical UF behavior 
of 4> w [x) at non-fouling conditions, namely its monotonic 
concave-shaped increase with increasing x caused by the 
steady-state advective-diffusive mass balance inside the 
CP layer. An asymptotic analysis based on Eqs. (28) and 


as 


<t>o 


1 + Ax 1/3 , 


(53) 


in the immediate vicinity of the fiber inlet (not re¬ 
solved on the scale of Fig. [5]), with A = 1.857 x 
T 2l 1/3 

v w/ 7^(^o) • Thus, (f> w (x ) shares the fractional x- 

dependence of the CP layer thickness S(x) for very small 
x/L . As expected, the initial growth of the surface 
concentration at the inlet is larger for larger TMP and 
larger membrane permeability L p (larger u 3 ) and smaller 
for larger characteristic shear rate 7 and bulk diffusion 
coefficient D(</> 0 ). Since the curves for (f> w (x ) are non¬ 
intersecting for different values of 7 , D{cj ) 0 ), and , these 
trends remain valid for non-snrall values of x/L. 

In some earlier UF works Ei ED], the concentration 
dependence of the collective diffusion coefficient, D(</>), 
for impermeable particles was accounted approximately, 
but the dispersion viscosity was treated as constant in¬ 
side the thin boundary layer, of value equal to the bulk 
(feed) viscosity r]((j ) 0 ). It is noticed from the blue dot¬ 
ted line in Fig. [ 5 J obtained for constant 77 = r]((f>o)i that 
at larger distance from the inlet, the growth of 4> w {x) 
with increasing x is underestimated in comparison to the 
surface concentration profile for variable viscosity, with 
the latter increasing for increasing concentration in ac¬ 
cord with Fig. [4} The reason for this is, as quantified in 
Eq. (22), that the axial particle transport lowering the 


CP layer formation is overestimated for a constant 77 in¬ 
side the layer taken equal to the minimal feed dispersion 
value. Commonly, 0o ^ 1 so that the feed viscosity and 
collective diffusion coefficient of electrically neutral parti¬ 
cle systems are practically equal to the fluid viscosity, 770, 
and single-particle diffusion coefficient, D 0 , respectively 
(see Figs. [3] and [4]). 

In most classical UF models, both D and 77 are de¬ 
scribed as being constant E1I27MIIII3 and set equal to 
the feed dispersion or infinite dilution values. The sur¬ 
face concentration profile (red dashed line in Fig. [5]) for 
constant D(tj) 0 ) and r)((f>o) is located in between that for 
concentration-variable D and 77, and for variable D((p) 
and constant 77 = 77(^0) (blue dotted line), respectively, 
for the reason that the transverse diffusion flux of parti¬ 
cles away from the membrane surface is underestimated 
when the minimal bulk dispersion value D(cf> 0 ) is used 
instead of a concentration-variable D(cf>) (see again Fig. 
[ 3 ]). The main part of Fig. [ 5 ] states that the precise 
modeling of the concentration dependence of D and 77 


(13) shows that the surface concentration grows initially dotted 


is key to a proper identification of operating conditions 
for which membrane fouling caused by the cake formation 
of jammed particles can be avoided. 

In the inset of Fig. [5j the similarity solution result 
for (f> w (x) with concentration-dependent D(</>) and 77 {(j>) 
(black solid line) is compared with two semi-analytic 
solutions for constant D and 77 of the transcendental 
Eq. ((36 ) (colored dashed-dotted lines). The red dashed- 
ine is the result of the transcendental equation 
for constant D{(f) w ) and ? 7 (</>o) (largest value of D , given 
at the membrane surface, and smallest possible viscosity 
value, given in the bulk). Since the extent of CP layer 
formation is underestimated in this way, a lower bound 
for the surface concentration profile with concentration- 
dependent transport coefficients is obtained. An up¬ 
per bound follows when the minimal diffusion coefficient 
D{(j>o) and the maximal viscosity are used as con¬ 

stant input values in Eq. (36) (blue dashed-dotted line). 


The strong increase of the blue dashed-dotted line at 
small x/L is a consequence of approximating 77 by its 
maximal membrane surface value y((f> w ) throughout the 
CP layer. 

Sufficiently distant from the inlet, the concentration 
profile for concentration-variable D and 77 is well de¬ 
scribed by the arithmetic average of the two bounding 
red and blue dashed-dotted lines (black dashed-dotted 
line). Since the numerical solution of Eq. ( 36 ) is fast 


and stable, the arithmetic average result is convenient 
for disclosing unwanted operating conditions where cake 
formation by particle excluded volume solidification takes 
place. 


B. Microgel permeability effects 

We proceed by discussing the influence of the particle 
permeability on the CP layer, cf>(x, y), and permeate flux, 
v w (x), profiles. The results shown in the following have 
been obtained from the numerical solution of the simi¬ 
larity solution scheme in Subsec. m As input, we have 
used the factorization approximation results for 77 (y, <(*), 
described by Eqs. (48)- © and pictured in Fig. [4 and 
the polynomial fits to the high-precision simulation data 
of D(x, 4 i), depicted in Fig. [3| 

Figure [ 6 ] shows how the membrane-surface concentra¬ 
tion profile, <j> w (x), is affected by the particle permeabil¬ 
ity (x = 20, 50, and 00 ), using two different particle 
radii a. We recall from Eq. (38) that y = 20, for exam¬ 


ple, amounts to a hydrodynamic penetration depth equal 
to the twentieth part of the particle radius. According 
to Fig. [ 6 j the membrane surface concentration decreases 
with increasing permeability (decreasing y). The rea¬ 
son for this is that for a given concentration D(x , <j>) in¬ 
creases and rj{4>) decreases with increasing permeability 
(see Figs. [3] and [4|. Both effects enlarge the number 
of particles leaving the CP layer, with <j> w being reduced 
accordingly. The diminution of 4> w (x) is most significant 
at the fiber outlet x = L. We also notice from Fig. [ 6 ] a 
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strongly enhanced surface concentration when the parti¬ 
cle radius is enlarged from 20 nm (red lines) to 30 nm 
(black lines). This can be mainly attributed to the par¬ 
ticle size dependence of the single-particle diffusion co¬ 
efficient, Do(x)> in Eq. (39), which is proportional to 
1/a for constant y. Recall here the [_D(y, <po)]~ 2 ^ 3 de¬ 
pendence of the coefficient A in Eq. (531, quantifying 
the surface concentration growth at the inlet. 

The transverse microgel concentration profile, <p(x,y), 
in the CP layer is shown in Fig. [7] as a function of the 
distance y from the membrane surface, for different re¬ 
duced axial positions x/L indicated in the figure. As 
described implicitly by Eq. (35), the transverse profile 
decays strictly monotonically from the membrane sur¬ 
face value, <f> w (x), at y = 0 towards the bulk concentra¬ 
tion value, <f> o, asymptotically reached for y >> 5(x). Due 
to the gradual buildup of the CP layer along the fiber, 
the transverse particle concentration profile is larger for 
larger axial distance x from the inlet. The effect of par¬ 
ticle permeability is to decrease the transverse concen¬ 
tration profile. This decrease is most pronounced at the 
membrane surface (dashed and solid lines). 

Results for (f> u ,(x) and the permeate velocity profile, 
v w (x), along the membrane surface are shown in Figs. 

? a) and (b), respectively, for different combinations of 
P and 7 as indicated, and y = 20 and oo (solid and 
dashed lines), respectively. The permeate velocity dis¬ 
played in Fig. [ 8 ] decreases monotonically with increasing 
x , due to the gradually sharpening CP layer structure: 
the particle concentration (f> w increases with the axial dis¬ 
tance from the inlet, thus enhancing the osmotic pressure 



FIG. 6 . Membrane surface concentration profile for perme¬ 
able particles for two particle radii a as indicated (black and 
red lines). Solid and dotted lines are similarity solution re¬ 
sults for inverse permeability parameters y = 20 and y = 50, 
respectively. Dashed lines represent results for impermeable 
hard spheres, y = oo. Other parameters: A P = 300 Pa, 
7 = 75 s“\ (f>o = 1.0 x 10 —3 , and L p = 5 x 10 -9 mPa - 1 s -1 . 



ys(L) 

FIG. 7. Transversal boundary layer concentration profile, 
at different axial positions x as indicated (differently 
colored lines). The transverse distance, y, from the membrane 
surface is scaled by the CP layer thickness at the fiber outlet, 
5(L). Solid lines: permeable particles with y = 20. Dashed 
lines: impermeable particles, y = oo. Other parameters as in 
Fig. 0 for a = 30 nm. 



x/L 


FIG. 8 . (a) Surface concentration profile, (p w (x ), and (b) per¬ 
meate velocity profile, v w (x), scaled by the pure solvent flux, 
v 9 , for different TMPs, A P, and shear rates, 7 , as indicated 
(differently colored lines). Solid lines: permeable particles 
with y = 20. Dashed lines: impermeable particles, y = 00 . 
Other parameters as in Fig. [5] 


at the membrane surface, which in turn, accord¬ 

ing to Eq. (13), lowers the values of v w (x) for given TMP. 
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FIG. 9. Fiber-length-averaged permeate velocity, (v w ), as 
function of the applied TMP for different shear rates, 7 , and 
feed particle concentrations, fa, as indicated (differently col¬ 
ored lines). Solid lines: permeable particles with x = 20. 
Dashed lines: impermeable particles, x = 00 ■ The dotted line 
applies for pure solvent where (v w ) = L p AP. Inset: fiber- 
length-averaged osmotic pressure, (II), scaled by the TMP. 
Other parameters as in Fig. [5] 


Since <f> w increases as x 1 ^ 3 in the immediate vicinity of the 
fiber inlet, and since II(0^) ex <p w + O(0£,) near the inlet 
where tf> w ~ fa <C 1, the permeate velocity decreases as 
a ; 1 / 3 very close to the fiber inlet. 

According to the black and blue lines in Fig. [ 8 ja) for 
A P = 5 kPa and 3 kPa, respectively, an enlargement of 
A P for fixed 7 causes an increase in <f> w due to the larger 
permeation drag experienced by the particles (larger v^ 
in Eq. ©)■ Quite interestingly, an increase of the TMP 
by 2 kPa is overcompensated in the ratio / AP in 

Eq. ( p~3| ) by the associated enlargement of n(^> xtJ ), caus¬ 
ing a decrease of the relative permeate velocity, v w /, 
for enlarged TMP as depicted in Fig. §b). While the 
relative permeate velocity is lowered, its absolute value, 
v w , is still enlarged for enlarged TMP. Different from 
this, if 7 , and thus the influx of the feed particles, in¬ 
creases for fixed A P (red and blue lines), the resulting 
enhanced axial drag gives rise to a smaller values of the 
membrane surface concentration and associated osmotic 
surface pressure and hence to a larger absolute permeate 
velocity. Notice here that 5(x) oc 7 -1 / 3 according to Eq. 
{±9). 

Figures [ 8 ja) and (b) illustrate, in addition, the influ¬ 
ence of the particle permeability: since (f> w decreases with 
increasing permeability, the permeate velocity, v w , at 
large distances from the inlet is for particles with y = 20 
significantly larger than for impermeable ones (solid ver¬ 
sus dashed lines). 

Consider next the fiber-length-averaged permeate ve¬ 


locity, 


{v w ) =-jr J v w (x)dx, 


(54) 


which according to Eq. (13) is related to the fiber-length- 
averaged osmotic surface pressure, (II), and the applied 
TMP by 


M 

AP 


= 1 - 


(v w ) 


(55) 


The average permeate velocity is more easily accessible 
experimentally than the local axial velocity distribution 
v w (x). Results for (v w ) as a function of AP are presented 
in Fig. [9] for different 7 and feed concentrations fa as in¬ 
dicated. For small values of the TMP up to AP « 3 kPa, 
the CP layer and the associated osmotic surface pressure 
are so weak that (v w ) ~ v^, = L p AP is valid. With fur¬ 
ther increasing TMP, the average permeate flux increases 
weaker than linearly in AP. This is explained in the in¬ 
set of Fig. [ 9 ] which shows the pressure ratio (II)/AP 
as a function of AP. The CP structure development 
along the membrane surface is such that (II) increases 
more strongly than linearly with increasing AP, giving 
rise to an average permeate velocity smaller than that 
for pure solvent feed, according to Eq. (131. Note fur¬ 


ther that (v w ) decreases with increasing fa since <j> w is 
increased. This effect becomes more significant with in¬ 
creasing AP. Moreover, consistent with the observations 
made in Fig. [ 8 ] for the axial permeate velocity distribu¬ 
tion, (v w ) increases with increasing 7 and with increasing 
particle permeability (dashed and solid lines). 

The absence of a plateau region for (v. w ) in Fig. [ 9 ] at 
large TMP indicates that the operation conditions for a 
limiting flux behavior at large TMP [43j have not been 
approached. Operating UF under limiting flux conditions 
reduces the efficiency of the separation process. The dis¬ 
cussion of (v w ) as a function of the applied TMP, such 
as the one in Fig. [9j is thus also of practical relevance. 

An important conclusion of the present subsection is 
that the particle permeability should be included in the 
UF modeling of non-ionic microgel dispersions. At large 
TMP, the permeate velocity (membrane surface con¬ 
centration) of permeable particles is significantly larger 
(smaller) than the one of impermeable particles, for oth¬ 
erwise identical system parameters. 


C. Efficiency and cost indicators 

To assess the performance of a filtration process from 
different technical-economical viewpoints at given oper¬ 
ating conditions, various efficiency and cost indicators 
can be introduced and evaluated. We discuss in the fol¬ 
lowing several indicators and their inter-relations. These 
indicators apply quite generally but are evaluated here in 
the context of the considered steady-state thin boundary 
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FIG. 10. Degree of Concentration Factor, a , as a function of 
applied TMP (a) and shear rate (b). Other parameters as in 
Fig. [5] 



FIG. 11. Productivity per Unit Membrane Area, 9, as a func¬ 
tion of applied TMP (a) and shear rate (b). Other parameters 
as in Fig. [5] 


layer cross-flow UF in a hollow cylindrical fiber mem¬ 
brane of length L and inner radius R. A more detailed 
discussion of the process indicators is given in the sup¬ 
plementary information. 

The first efficiency indicator in the present list of in¬ 
dicators is the Degree of Concentration Factor , a > 1, 
defined as the ratio of the particle volume fraction, iff, 
in the finally obtained dispersion (retentate) to the feed 
dispersion volume fraction, 


h_ = 1 

</>o 1 - P ‘ 


(56) 


Here, /? < 1 is the related so-called Solvent Recovery 
indicator which, by mass and dispersion volume conser¬ 
vation, is equal to the fraction of initial dispersion volume 
recovered in the permeate compartment as pure solvent. 
For fully developed homogeneous parabolic Poiseuille 
feed flow described by Eqs. <© and ( p~7] ) , the second 
indicator is obtained in the supplementary information 
as 


P = 


4 L ( v w ) 

R O'm. 


(57) 


where u m = 7 P /2 is the feed inflow velocity at the center 
of the fiber. 


According to Eq. (57), /? is proportional to ( v w ), thus 
having an explicit 1/7 dependence inherited from u m . 
Considering the TMP and characteristic shear rate de¬ 
pendence of ( v w ) displayed in Fig. [9j it follows that /3 in¬ 
creases with increasing AP, first linearly for small TMP 
values and subsequently sub-linearly for larger pressure 
values. In contrast, /3 decreases monotonically with in¬ 
creasing 7 . Since a increases with increasing /3 as quan¬ 
tified by Eq. (56), the Degree of Concentration Factor 


a increases with increasing AP and decreases with in¬ 
creasing 7 , as expected. Figs. [l0](a) and (b) show quan¬ 
titative results for a as a function of AP and 7 , respec¬ 
tively, obtained from the similarity scheme solution for 
concentration-dependent transport properties. 

An important indicator for the material cost aspect of 
a baromembrane separation process is the Productivity 
per Unit Membrane Area indicator, 9, defined as the dis¬ 
persion volume flux at the outlet (retentate flux) divided 
by the membrane area. It can be expressed in terms of 
a and (v w ) as 


a 2 (v w ) 


a — 1 


(58) 


We refer to the supplementary information for the 
derivation of this equation. Note that 9 has the phys¬ 
ical dimension of a velocity. An alternative definition of 
the process productivity using the permeate rather than 
the retentate as the final product is discussed in Ref. 69]. 

As shown in Fig. |ll[ a), 9 increases monotonically with 
increasing AP, with a more pronounced upturn for the 
smaller considered shear rate. Quite interestingly, the 
behavior of 9 as a function of 7 is non-monotonic, with 
a minimum value at a shear rate that increases with in¬ 
creasing TMP (Fig. &))• Increasing the shear rate is 
associated with two opposite trends. The first one is a 
reduction in the CP layer that results in a decrease of 
the osmotic pressure opposing the TMP. This leads to 
a diminution of 9. The second trend is an increase in 
the feed flow that leads to an increase in the outlet flow, 
with a correspondingly enlarged indicator 9. At suffi¬ 
ciently low shear rates, the osmotic pressure noticeably 
compensates the TMP. Therefore, the first noted trend 
predominates thus lowering the CP layer profile accom¬ 
panied by a stronger permeate flux and a lower outlet 
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flow. When the shear rate is sufficiently high, the CP 
layer profile is reduced to such an extent that its further 
reduction does not substantially affect the permeate flux 
any more. Hence, the second noted trend leads to an 
increase in 9 with increasing 7 . 

We now discuss two energy cost indicators. The first 
one, termed Specific Energy Consumption indicator, w, is 
defined as the energy consumed for producing a unit vol¬ 
ume of final product. In earlier related works EDHISj, the 
final product has been identified with the permeate (for 
example, with desalinated water in case of desalination 
studies) whereas the final product is identified here with 
the retaining dispersion of concentration iff. Moreover, 
in earlier works the CP layer has been either completely 
disregarded or treated in a simplified way using a semi- 
empirical parameter m ■ Considering that basically the 
whole externally supplied power is consumed for press¬ 
ing the solvent through the membrane, it is shown in the 
supplementary information that u> can be expressed as 

w = (a-l)AP, (59) 


so that its shear rate dependence is inherited from that 
of a. Consequently, uj decreases with increasing 7 . More¬ 
over, it strongly increases with increasing AP, quadrati- 
cally in the pressure for small applied TMP or high shear 
rates, and more stronger than linearly but weaker than 
quadratically for larger TMP values. 

While the indicators a, 9 , and oj discussed thus far fa¬ 
vor large values of the applied TMP for optimizing the 
one-stage cross-flow UF process, this conflicts with the 
second energy cost indicator considered here, namely the 
indicator e < 1 of Specific Energy Efficiency. This indica¬ 
tor is defined m as the ratio of the thermodynamically 
necessary minimal work per unit volume of final product, 


Wmin, to the energy per volume, w, in Eq. (591, consumed 
during the steady-state filtration process, 


oufo 

(a - 1) 


cm s- » 


Explicitly, 


^min 


r v f r<t>f 

/ n dv = <t> } l n(</>) 

JVq j 4> 0 


1 

VfJvo 


d<f> 


(61) 


is the minimal work per unit volume of retentate required 
in a adiabatically slow reversible thermodynamic process 
of isothermally compressing the dispersion from its initial 
volume Vo to the final volume Vf with the help of a fully 
particle-retentive non-adsorbing membrane employed as 
a plunger. We have used the relation V = ( 47 r/ 3 )a 3 ./V/</>, 
with N the total number of mobile particles that remains 
constant during the compression. Note that the minimal 
work required in the adiabatically slow reverse-osmosis 
compression is equal to the change in free energy of the 
system minus the energy required for lowering the num¬ 
ber of solvent molecules in the system during the com¬ 
pression. The integral in Eq. (60) can be performed 



FIG. 12. Specific Energy Efficiency, e, as a function of applied 
TMP (a) and shear rate (b). Other parameters as in Fig. [ 5 ] 


analytically, using the CS Eq. (371 as input for the os¬ 
motic pressure n((/>). The resulting expression for e is 
given in the supplementary information. 

The behavior of the Specific Energy Efficiency as a 
function of TMP and characteristic shear rate is shown 
in Figs. |T2j(a) and (b), respectively. As expected, e de¬ 
creases with increasing 7 , but different from the earlier 
discussed indicators, e decreases with increasing AP. It 
should be noticed here that while the Degree of Con¬ 
centration Factor a increases with increasing TMP, the 
factor a /(1 —a) in Eq. (60) decreases. We finally point to 
the small values of the Specific Energy Efficiency in Fig. 


12 This is due basically to the comparatively large par¬ 


ticle hard-core radius a = 10 nm for which the osmotic 
pressure buildup along the membrane surface is small as 
described by the size dependence n oc 1 /a 3 . 


V. CONCLUSIONS 

We have presented the first comprehensive theoretical 
study of particle permeability effects in the UF of non¬ 
ionic microgel dispersions modeled as solvent-permeable 
Brownian spheres. Based on a boundary layer analysis, 
we have provided a systematic study of inside-out cross- 
flow UF. The CP layer and permeate flux profiles have 
been calculated using accurate analytic expressions for 
the concentration-dependent collective diffusion coeffi¬ 
cient, D(f>), and zero-frequency viscosity, of disper¬ 
sions of permeable particles. The employed expressions 
for these transport coefficients are well tested against 
computer simulation results and experiments j40l . They 
should be used instead of the substantially less accu¬ 
rate cell model results and the semi-empirical Krieger- 
Dougherty and Eilers expressions employed in earlier the¬ 
oretical UF studies, which have dealt mostly with imper- 
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meable colloidal particles and proteins. 

We have shown that the concentration dependence 
both of D and 77 significantly influences the CP layer, in 
particular right at the membrane surface, with a corre¬ 
spondingly strong influence on the osmotic pressure pro¬ 
file and therefore on the permeate velocity. The usage 
of constant values for D and 77, if taken in particular for 
the common situation of low feed volume concentration 
(f> 0 , results in the underestimation of the CP layer struc¬ 
ture and hence in the overestimation of the permeate 
flux. The applicability of the transcendental Eq. (361 for 
the membrane surface concentration profile, < p w (x ), with 
constant transport coefficients has been discussed (inset 
in Fig. [ 5 ]). It can be used to generate an arithmetic aver¬ 
age concentration profile in good agreement with calcu¬ 
lations of 4> w {x ) for concentration-dependent D and 77 at 
non-small distances x from the fiber inlet. Since the nu¬ 
merical solution of the transcendental equation is quick 
and stable, the arithmetic-average profile can be prof¬ 
itably used, for example, in determining the axial range 
extending downstream from the inlet wherein the mem¬ 
brane is devoid of a stagnant cake layer [nmsiET]. 


We have shown that the effect of non-zero microgel 
permeability is to reduce the particle concentration at 
the membrane surface with a consequently enlarged per¬ 
meate velocity. It is worth to emphasize again that our 
study of permeability effects is relevant not only for non¬ 
ionic microgels, modeled here hydrodynamically as uni¬ 
formly permeable spheres, but also for other hydrody¬ 
namically structured particles such as core-shell parti¬ 
cles with a dry rigid core and an outer permeable poly¬ 
mer brush layer. As pointed out earlier and discussed in 
Ref. [30]) the reason for this is that even particles with 
a complicated internal structure are hydrodynamically 
fully characterized, for quite general conditions, by their 
hydrodynamic radius. 

Various efficiency and cost indicators assessing 
the steady-state UF process from different technical- 
economical viewpoints have been discussed in depen¬ 
dence of the process parameters and with full account 
of the CP layer structure. A detailed description of the 
indicators and their inter-relations is given in the supple¬ 
mentary information. Different indicators favor partially 
conflicting optimization strategies. For example, while an 
enlarged TMP leads to a larger Degree of Concentration 
factor, a, and Productivity per Unit Membrane Area, 9 , 
the price to pay is a reduced Specific Energy Efficiency, 
e. The discussed indicators are of value going beyond 
the specific permeate-particles cross-flow UF process dis¬ 
cussed in the paper. They can be used, for example, 
to evaluate different optimization strategies employing 
multistage processes (TDHza, re-use of free energy accu¬ 
mulated in a filtration process [73HZ5|, external electric 
fields [76] , local production of flow instabilities [78], and 
the generation of high shear rates by moving parts m- 
We expect our theoretical work to be useful in future 
experimental explorations of microgel cross-flow UF. 

The present paper was focused on identifying and ana¬ 


lyzing hydrodynamic effects on UF arising from the par¬ 
ticle permeability and concentration-permeability depen¬ 
dencies of D and 77, with fouling effects disregarded. The 
dynamics of cake layer formation by permeable parti¬ 
cles can be described using cross-sectional averaged forms 
of the time-dependent advection-diffusion and stationary 
momentum balance equations dn na m eq], with the 
local hydraulic resistance of the cake layer added to the 
clear membrane resistance R m . 

In dealing with UF of microgel dispersions, also par¬ 
ticle polydispersity and deformability play a role. Size 
polydispersity is expected to enlarge the particle concen¬ 
tration in the cake and CP layers. Moreover, there should 
be an enhanced tendency of membrane pore clogging by 
the small particles of the size distribution. There is also a 
tendency of size fractionation triggered by the transver¬ 
sal flow which, however, is counteracted by the strong 
Brownian motion in UF. The theoretical treatment of 
polydispersity in UF is a difficult task since it requires the 
additional consideration of collective cross-diffusion coef¬ 
ficients relating the diffusion flux of one component to the 
local concentration gradient of another one [321 1ST! 82] , 
These coefficients are affected, in particular, by His caus¬ 
ing a hydrodynamic dragging-along of small particles by 
larger ones. 

For ionic microgels in polar solvents, the effective par¬ 
ticle and membrane surface charges are additional im¬ 
portant parameters j65] [53] [51] . As shown in Ref. [T] , 
while the permeability of ionic microgels is hydrodynam¬ 
ically less influential than that of uncharged ones, due 
to the longer-ranged electrostatic repulsion, the ionic mi¬ 
crogels tend to shrink substantially with increasing con¬ 
centration in the swollen-state temperature range. More¬ 
over, the softness (single-particle compressibility) of ionic 
microgels influences the width of the fluid-solid coex¬ 
istence region, rendering it smaller for stiffer microgels 
(85l . Thus, the size shrinkage with increasing concentra¬ 
tion and the softness of ionic microgels must be included 
in a realistic modeling of the CP and cake layers. Work 
by the present authors on filtration modeling of polydis- 
perse and charge-stabilized rigid colloidal particles and 
soft ionic microgels is in progress. Furthermore, we are 
currently involved in an experimental-theoretical study 
of cross-flow filtration of microgel dispersions [86]. The 
present theoretical results in particular for the concen¬ 
tration factor and energy efficiency indicators are used 
to analyze the performance and to design the filtration 
experiments. The outcome of this ongoing study will be 
communicated in a future publication. 
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S-l 


SUPPLEMENTARY INFORMATION: 

PROCESS INDICATORS 

We discuss here in more detail the process indicators 
introduced in Subsec. IIV Cl These are functions of the 
adjustable operation parameters and here in particular 
of the volume fraction, (f> o, of colloidal particles in the 
homogeneous feed dispersion. For inside-out cross-flow 


UF in a hollow cylindrical fiber membrane (see Fig. SI) 
the indicators are obtained from the calculated spatial 
distributions of the local particle volume fraction, cf>{r , z), 
and the dispersion velocity, 


v(r,z) = u r (r,z)e r +u z (r,z)e z 


(S.l) 


expressed in cylindrical coordinates, (r, z), with the ra¬ 
dial velocity component u r along the unit vector e r and 
the axial velocity component u z along the unit vector e z . 
In this global coordinate frame, the velocity components 
at the inner membrane surface, r = R, are u r (R,z) = 
u(x,y = 0) and u z (R,z) = v(x,y = 0) = —v w (x) for 
x = z, where u(x, y) and v(x, y) are the longitudinal and 
axial velocity components measured, respectively, in the 
membrane surface anchored local coordinate frame used 
in our CP layer analysis (see Eq. ([?]) and Fig. [2]). 


For a steady-state UF process where V • J = 0 is valid 
according to Eq. |2]) and for a fully particle retentive 
membrane where Eq. (18) applies to, the longitudi¬ 
nal particle fluxes through different fiber cross-sections 
should be equal. In particular, 

f dSc/>(r, L)u z (r, L) = (j> 0 f dSu z (r,0) (S.3) 

Js Js 

where we have ignored once again the small longitudinal 
diffusion flux contribution. The right-hand-side of Eq. 


(S.3) is the inlet particle flux of the feed dispersion of 


uniform volume fraction </>(r, 0) = </> 0 . By combining Eqs. 
(S.2) and (S.3), the expression 


<Pf _ Is dSu z (r, 0) 




(f> 0 f s dSu z (r,L) 


(S.4) 


for the Degree of Concentration factor a is obtained, 
equal to the ratio of inlet to outlet dispersion volume 
flux. 

Another expression for a, particularly suitable for the 
presented boundary layer analysis, is obtained from dis¬ 
persion volume conservation (V • v = 0), 


f dSu z (r } 0) = [ dSu z (r,L) 

Js Js 
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FIG. SI. Global cylindrical coordinate system ( r,z ). The z- 
coordinate line extends along the axis of the hollow cylindrical 
fiber membrane of length L and inner radius R. 


1. Final product efficiency 

The final product is identified in the present work 
with the retentate dispersion of mean volume concen¬ 
tration <j)f, flowing through the outlet cross-section at 
z = L into a final reservoir. The particle advective flux 
through the outlet cross-section is f s dS c/)(r, z)u z (r, L) 

where f s dS = 2n rdr. Hence, <fif is given by the ratio 
of the longitudinal particle advective to the longitudinal 
dispersion volume flux, 


cj)f = 


J s dSL)u z (r, L) 
J s dS u z (r,L) 


(S.2) 


where we have neglected the very small longitudinal dif¬ 
fusion flux contribution. 


+ 27tI? / dzu r (R,z), (S.5) 

Jo 

stating that the volume inflow through the inlet cross- 
section is equal to the sum of the volume outflows 
through the outlet cross-section and the cylindrical mem¬ 
brane of area 2nRL. The combination of Eqs. (S.41 and 


(S.5) yields 


“=rb- 

with the Solvent Recovery indicator, 

Jo dzu r (R , z) 


P = 


2ttRL L 


J s dSu z (r, 0) 


(S.6) 


(S.7) 


equal to the fraction of initial dispersion volume recov¬ 
ered in the permeate compartment. Using u r (R,z) = 
—v w (x) together with Eq. (54) for (v w ) and the cylindri¬ 
cal Pouiseuille flow profile for u z (r, 0), the Eq. (57) for (3 
in the main text is obtained. 

The logical way of calculating the final product effi¬ 
ciency indicators is to determine first /? using Eq. (S.7), 
with a and <determined subsequently using Eq. (S.6). 


2. Productivity per unit membrane area 

The indicator of Productivity per Unit Membrane 
Area, 0, is defined here as the retentate flux divided by 
the membrane area, 


6 = 


f s dS u z (r, L) 


2-kRL 


(S.8) 
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giving this indicator the dimension of a velocity. Using 
Eqs. (S.4)-(S.6), 9 can be related to (3 by 


9 = 


J s dSu z (r, 0) _ f Q dzu r (R,z) 
(1 - /3)2ttRL ~ /3(1-P)L 


(S.9) 


The integral in the second equality is proportional to 
the fiber-length-averaged permeate velocity. Using Eqs. 
(54) and (56) in the previous expression, Eq. (581 in the 


main text is obtained. 


3. Energy cost 


We discuss first the indicator of Specific Energy Con¬ 
sumption, w, defined as the energy consumed to produce 
a unit volume of final product which, in our case, is the 
retentate dispersion of volume fraction <pf. Thus, for a 
steady-state process, w is equal to the ratio of consumed 
power to the outlet volume per unit time. Assuming that 
basically the whole external power is spent in pressing the 
solvent through the membrane, this implies 


2-kR fg dz AP(z)u r (R, z) 
J s dSu z (r,L) 


(S.10) 


where A P(z) is the local transmembrane pressure differ¬ 
ence which in general depends on the axial distance z 
from the inlet. In our calculations, the TMP has been 
taken as constant. Note that the physical dimension of 
oj is energy per volume. 

In general, one can expect that variations of A P(z) 
along the hollow fiber membrane are much smaller than 
the mean pressure value, denoted here as AP, so that 


A P(z) approximated by A P can be shifted out of the 
integral in Eq. (S.lOl. Consequently, by using in addi¬ 
tion Eqs. (S.41 and (S.61, Eq. (591 in the main text is 
obtained wherein lo is expressed by the product of the 
fiber-length-averaged TMP A P and (a — 1). 

As discussed in Subsec. |IV C[ the Specific Energy Effi¬ 
ciency indicator, e, is defined as the ratio of the thermo¬ 
dynamically necessary minimal reverse-osmosis compres¬ 
sion work, w m i n , required to produce a unity of retentate 
dispersion to the Specific Energy Consumption ui. Ac¬ 
cording to Eq. (601, e can be expressed in terms of (j> 0 , 
AP, and osmotic pressure 11(2) along the membrane sur¬ 
face as 


a <j>o 


r<h 


{a - 1) A P 


d(f> 


vo 


m 


(S.ll) 


Substituting II(</>) as described by th e CS Eq. (37) 
for colloidal hard spheres into Eq. (S.ll), we obtain the 
analytic expression 


3 ksT acf> o 

47 ra 3 (a — 1) A P 


In (a) + 


3 — 2 cj)f 3 — 2<f>o 


(i -h) 


(1 - </>o) 
(S.12) 


for e in terms of the input parameters 4> o, AP, and hard¬ 
core radius a, and also in terms of a = <j>f/(j )o which 
must be determined p rior t o e. The factor to the left 
of the bracket in Eq. (S.12) can be recast into the sug¬ 
gestive form ( k B Tnb/AP ) (a /{a — 1)) invoking the ratio 
of the van’t Hoff pressure ksTrib and the TMP, where 
rib is the number density of the injected feed dispersion. 
Note further the explicit 1/a 3 dependence of e which ren¬ 
ders this indicator very small for larger colloidal particles 
where the osmotic pressure buildup along the membrane 
is small. 



































